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We calculate how strongly one can constrain the alternative theories of gravity with deci- 
Hz gravitational wave interferometers such as DECIGO and BBO. Here we discuss Brans- 
Dicke theory and massive graviton theories as typical examples. We consider the inspiral 
of compact binaries composed of a neutron star (NS) and an intermediate mass black hole 
' (IMBH) for Brans-Dicke (BD) theory and those composed of a super massive black hole 

fS) , (SMBH) and a black hole (SMBH) for massive graviton theories. Using the restricted 2PN 

waveforms including spin effects and taking the spin precession into account, we perform the 
Monte Carlo simulations of 10* binaries to estimate the determination accuracy of binary 
parameters including the Brans-Dicke parameter ojbd and the graviton Compton length Ag. 
Assuming a (1.4, 1Q)Mq NS/BH binary of SNR=\/200, the constraint on ojbd is obtained 
' as ojbd > 2.32 x 10^, which is 300 times stronger than the estimated constraint from LISA 

observation. Furthermore, we find that, due to the expected large merger rate of NS/BH 
binaries of 0(10*) yr~^, a statistical analysis yields ujbd > 3.77 x 10*, which is 4 orders 
^ I of magnitude stronger than the current strongest bound obtained from the solar system 

>j» ' experiment. For massive graviton theories, assuming a (10", 1O'^)M0 BH/BH binary at 3Gpc, 

}_( ' one can put a constraint Ag > 3.35 x lO^^cm, on average. This is three orders of magnitude 

{DJO, stronger than the one obtained from the solar system experiment. From these results, it is 

understood that DECIGO /BBO is a very powerful tool for constraining alternative theories 
of gravity, too. 
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Many challenges have been made to modify gravitational theory from general 
relativity in order to explain the current acceleration of the universe.l^ In this letter, 
CN ' we focus on two simple possibilities, Brans-Dicke theor}® and massive graviton the- 

, ories.'^'"^' Brans-Dicke theory is a scalar-tensor theor}^^ of the simplest type. This 

OO \ theory is parameterised by the so-called Brans-Dicke parameter ljbd and approaches 

^ ■ to general relativity in the limit wbd — ^ co. The current strongest bound on this pa- 

Q , rameter, cjbd > 4 x 10"^, is obtained by the solar system experiment using the Saturn 

probe satellite Cassini.-' It measured the Shapiro time delay, which corresponds to 
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measuring the spatial metric deviation from general relativity. Another constraint 
on wbd has been obtained from orbital period decay rate of a binary PSR J1141- 
d ' 6545, which consists of a neutron star and a white dwarf.'^' Since there exists scalar 

dipole radiation in Brans-Dicke theory, the orbital evolution changes from the one 
in general relativity, whose information is implemented in the orbital period decay 
rate. Bhat et al.-- found that constraint becomes Oq = 2a;gp+3 < 2.1 x 10^^ for 
Brans-Dicke theory, which corresponds to wbd > 2.38 x 10^. Although this bound is 
2 times weaker than the Cassini one, it has distinct meaning since it gives a direct 
constraint on scalar dipole radiation. 

On the other hand, massive graviton theories by definition introduce a finite 
mass rUg to the graviton. Verification of Kepler's third law in the solar system 
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experiment puts the lower bound on the graviton Compton length, = h/nigC, as 
Xg > 2.8 X W^'^cm^ 

Estimate for the constraint on alternative theories of gravity using gravitational 
waves from inspiral compact binaries has been studied in several papers !^~^ Re- 
cently, Berti et a/A^- calculated the constraint on wbd and with LISA,!^''^^ 
using restricted 2PN waveforms and performed Monte Carlo simulations. Following 
their analysis, we improved their calculation by including the spin-spin coupling a, 
small orbital eccentricity, and spin precession effects.'^ Under the assumption of 
the so-called simple precessior^^^^ and when we restrict our calculation for cir- 
cular binary orbits, we obtained the constraints on u^d and A^ as wbd > 6944 
using a (1.4, 1000)Mq NS/BH binary of SNR=V200 and A^ > 4.86 x lO^icm us- 
ing a (1O^,1O6)M0 BH/BH binary at 3Gpc. When we include eccentricities, we 
found that the constraint on wbd is unaffected as long as we include prior infor- 
mation whilst the one on A^ becomes A^ > 3.10 x lO^^cm. Therefore we can say 
that the effects of eccentricities are not so strong for both cases. At the same time, 
Stavridis and WilP^ estimated the constraint on A^ for a circular BH/BH binary 
including both spins of binary objects and taking spin precession into account. For 
a BH/BH binary of (10*^, 1O'^)M0 at 3Gpc, they obtained Ag > 5 x lO^^cm when 
the spin-spin precession effect is taken into account and A^ > 4 x 10^^ cm when it 
is not taken into account. To compare our results with their ones, we estimated 
the constraint on Xg with a (10^, 1.1 x W^)Mq circular BH/BH binary under simple 
precession, in which the spin-spin precession effect is neglected, and obtained the 
constraint A^ > 3.7 x lO^^cm.'^^ Although we cannot directly compare these two 
results, it seems that our results are consistent with the ones obtained by Stavridis 
and Will.'^ There is also a recent work done by Yunes and Pretoriu^^ in which 
they proposed a new framework, the parametrised post-Einsteinian framework, to 
perform gravitational wave tests of alternative theories of gravity. 

Following our previous paper,!^ we estimate the possible constraint on wbd and 
Xg obtained by detecting gravitational waves from the inspiral of precessing com- 
pact binaries using deci-Hz space-borne gravitational wave interferometers such as 
DECIGd^^^S g^j^fj Bgo.l^''^ These detectors are most sensitive in the frequency 
band between O.lHz and IHz, and the noise levels are about four orders of magnitude 
lower than that of LISA. These detectors have a huge number of compact binaries 
as promising sources. (Because of that, high precision cosmology using them is also 
expected.'^Sl') In this letter we perform the analysis assuming DECIGO noise curve 
but almost the same results will apply to BBO, too. Here, we restrict our attention 
to circular binaries since the lower bounds on ljbd and A^ are not much affected by 
the inclusion of eccentricity into binary parameters. (The upper bound on lobd, if 
detected, can be affected by including eccentricity, though.) Since the detection rate 
of NS/BH mergers is expected to be 0(10"^) for DECIGO/BBO, we add a statistical 
analysis, which improves the constraint on the deviation from general relativity. Our 
results are only approximate estimates in that we do not take the errors coming from 
the use of approximate waveforms into accounlP^ and also due to the limitation of 
the validity of the Fisher analysis. 

First we review the basic plan of DECIGO (DECi-hertz Interferometer Grav- 
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itational wave Observatory) and show the noise curve. DECIGO is a planned space 
gravitational wave antenna mission .'^''^S it consists of four constellations of three 
drag free satellites which keep the triangular shape throughout the flight. These three 
satellites form Fabry-Perot interferometers with separation 1000km. DECIGO can 
detect gravitational waves from both astrophysical and cosmological sources mainly 
between 0.1-lOHz. Among them, compact binaries are the promising sources, al- 
though the primary goal of the mission is to detect primordial gravitational wave 
background at the level of i^cw = 10~^^. For LISA frequency band, such small sig- 
nals are completely obscured by the foregrounds generated by white-dwarf binaries. 
Since this foreground noise has cut off frequency around 0.2IIz,ED DECIGO has a 
much better chance to detect the primordial gravitational wave background. 
The instrumental noise spectral density for DECIGO is given bj® 



ST\f) = 5.3 X 10-^s 



2 2.3 X lO"'^ 2.6 X 10" 



where x = f / fp with fp = 7.36IIz. The three terms on the right hand side represent 
the shot noise, the radiation pressure noise and the acceleration noise, respectively. 

Besides instrumental noise, there are three main foreground confusion noises: 
galactic white dwarf binaries 5'^^'(/),'33i extra-galactic white dwarf binaries S'f^~^'^^ (f)^ 
and neutron star binaries S^^{f)^^ The noise spectral densities for the first two 
confusion noises are given in Ref. [T8|) . Following Ref. [MT) , the noise spectral density 
for NS binaries is estimated as S^^{f) = 1.3 x lO'^^^ (//IHz)"^/^ , where we have 
set the cosmological parameters to = 0.3, l7/i = 0.7 and Hq = 72km/s/Mpc. 

Then, the total noise spectral density for DECIGO becomes 



sr\f) 



^ (1) 
where the factor 7^(/) = exp{-2 (//O.OSHz)^} represents the cutoff of the WD/WD 
binary confusion noise. We put the factor 0.01 in front of which repre- 

sents the fraction of gravitational waves that cannot be removed after foreground 
subtraction. Namely, we assume that 99% of the gravitational waves from neutron 
star binaries can be identified and removed. Nf is the number density of galactic 
white dwarf binaries per unit frequency given in Ref. llSp . k ~ 4.5 is the average 
number of frequency bins that are lost when each galactic binary is fitted out and 
T represents the observation time which we fix as lyr. The noise curve ([T]) is shown 
in Fig. [1] as a thick solid curve. We introduce low and high cut-off frequencies at 
/low = 10~^Hz and /high = lOOHz, respectively. Although it has not been estimated 
rigorously, it might be possible to extend the noise curve on lower frequency side 
down to lO-'^-lO^^Hz.SS This bound comes from the limitation of controlling the 
mirror positions. This extension is shown as a thick dashed curve in Fig. [TJ We 
also show the noise curve of LISA as a thin solid curve for comparison. Thin dotted 
line and thick dotted line each, respectively, represent the amplitude of gravitational 
waves from a (1.4, 1OOO)M0 and a (1.4, 1O)M0 NS/BH binary of SNR p = 10 with 
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lyr observation before coalescence. Each dot labeled "lyr" represents the frequency 
at lyr before coalescence. 
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Fig. 1. The noise spectral density for DECIGO (thick solid curve) and LISA (thin solid curve). It 
might be possible to extend the DECIGO noise curve down to 10~^IIz^^' (thick dashed curve). 
We also show the amplitudes of gravitational waves from a (1.4, 1000)AfQ (thin dotted line) and 
a (1.4, lO)A/0 (thick dotted line) NS/BH binary with each SNR set to p = 10 for lyr observation 
before coalescence. Each dot labeled "lyr" represents the frequency at lyr before coalescence. 

We used the matched filtering analysis to estimate the determination accuracy 
of the binary parameters. There are 14 parameters in total: the chirp mass InA^, 
the dimensionless mass ratio r/; the coalescence time to the coalescence phase ^c! the 
distance to the source D\ the spin-orbit and spin-spin coupling coefficients, /3 and cr; 
the directional cosine between the orbital and spin angular momenta k, the precession 
angle parameter Qc; the direction of the source (^Sj'/'s); the initial direction of the 
total angular momentum <?!)j); finally, or A~i, depending on which theory we 
are aiming to constrain. We calculate the inverse of the Fisher matrix numerically, 
assuming stationary and Gaussian noise. Integration is performed from /in to /fin, 
where /in = max{/iow, /lyr} and /fin = min{/high, /isco}, respectively, hyr is the 
frequency at lyr before coalescence and /isco is the one at the innermost stable 
circular orbit (ISCO). First we performed the pattern-averaged estimate of the errors 
in determination of the parameters for binaries with various masses, in which we have 
averaged over the directions of the source and the orbital angular momentum, and the 
spins of the binary constituents are assumed to be zero. We set SNR=10 for Brans- 
Dicke theory and Dl = 3 Gpc for massive graviton theory. We also performed the 
Monte Carlo simulations for 10^ binaries, distributing {6s,4's,(^J,4'J^ ^) Q'c) randomly, 
both with and without the precession effect. We set the fiducial values to tc = 4>c = 
^BD ~ ~ ^ with (mNs^n-BH) = (1.4, 1O)M0 of SNR=10 (for pattern-averaged 
estimate) or SNR=-v/200 (for Monte Carlo simulations) for Brans-Dicke theory, and 
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Table I. The results of error estimation with DECIGO in Brans-Dicke theory for various mass 
NS/BH binaries. We performed pattern-averaged estimates, using only one detector and fixing 
the SNR to 10. 



masses 


/in 


/fln 




AhiM 


Alnr] 


Ap 


Aa 




(Hz) 


(Hz) 


(10') 


{10-') 


(%) 






(1.4, W)Mq 


0.118 


100 


1.342 


0.978 


2.78 


0.190 


2.18 


(1.4,5O)M0 


0.0776 


85.6 


0.2662 


2.34 


2.64 


0.106 


1.09 


(1.4, 100)Mq 


0.0651 


43.36 


0.1899 


2.34 


1.87 


0.0485 


0.563 


(1.4,400)Mq 


0.0460 


10.95 


0.04244 


4.96 


1.85 


0.0133 


0.250 



Table II. Comparison of the constraints on cjbd and other physical quantities for a (1.4, 1O)M0 
binary with DECIGO and a (1.4, 1OOO)M0 binary with LISA. We performed pattern-averaged 
analyses and SNRs are fixed to 10 for both cases. 



masses and detector 




uncor 


iVcw 




./in 


/fln 




(10«) 


(10«) 


(10«) 




(Hz) 


(Hz) 


(1.4, 10)Mq, DECIGO 


1.34 


332 


5.9 


0.027 


0.118 


100 


(1.4, 1OOO)M0, LISA 


0.00821 


21.6 


1.8 


0.083 


0.0366 


1.00 



("^bHj^ibh) = (1O^,1O^)M0 at Di = 3Gpc for massive graviton theories. For the 
analysis without the spin precession effect, the fiducial values of /3 and a are set to 
0. When we take it into account, we set the dimensionless spin parameter x = S/rn? 
to and 0.5 for the lighter and heavier bodies of binaries, respectively, where S is 
the magnitude of the spin angular momentum. We include prior information on (3 
and (T when calculating Fisher matrices. See Ref. jlSj) for more details. 

In Table HI we show the pattern-averaged results of binary parameter estima- 
tion errors for Brans-Dicke theory with (1.4, 10)Mo, (1.4, 50)Mo, (1.4, 1OO)M0 and 
(1.4,4OO)M0 NS/BH binaries of SNR=10. (It seems that SNRs of 0(10) are too 
small in performing the Fisher analysis.'® However, constraints for higher SNR 
binaries are obtained by just scaling in proportional to SNR. ) It can be seen that 
smaller mass binaries give stronger constraints on wbd- This can be understood 
as follows. The velocities of binaries at 1 yr before coalescences are slower for 
smaller mass binaries. Since Brans-Dicke theory gives dipole correction to binary 
gravitational waves, this correction is -IPN order. Therefore, when we fix SNRs, 
this contribution is larger for slower binaries, which makes the constraints stronger. 
Comparing these results with the ones in Ref. I18|) . we can say that DECIGO has 
better ability in constraining wbd compared to LISA. To be more exploit, let us com- 
pare the constraint from (1.4, 1O)M0 with DECIGO and the one from (1.4, 1OOO)M0 
with LISA. This mass parameter is an optimised choice for each detector to constrain 
Brans-Dicke theory. First, we compare the uncorrelated constraint Wg^"'^ which is 
calculated directly from the Fisher matrix (not the inverse of it) and represents the 
possible constraint when the degeneracies between wbd and other parameters have 
been solved completely.'^ Table IIIl shows that the former constraint is stronger than 
the latter by more than 1 order of magnitude. There are 2 reasons for this, (1) the 
number of gravitational wave cycles A'^qw are greater and (2) the velocity at 1 yr 
before coalescence viyr is slower. From Table [III we can see that Agw for the former 
case is about 3 times greater compared to the latter case. On the other hand, f lyr 
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Table III. The results of error estimation in massive graviton theory for BH/BH binaries at 3 Gpc 
with various masses. We performed pattern-averaged estimates using only one detector. 



masses 


/in 

(mHz) 


/fin 

(mHz) 


SNR 


A. 
(lO^^cm) 


AlnM 

(%) 


Aln-q 


Ap 


Aa 


(10^10^)MQ 


1.0 


2.20 


1338 


1.014 


14.0 


2.46 


9.40 


2.50 


(1O^1O^)M0 
(1O^1O5)M0 


1.0 


4.00 


2044 


1.270 


1.19 


1.69 


9.10 


2.44 


1.0 


22.0 


4909 


1.133 


0.0286 


0.930 


7.00 


1.69 


(lo^lo*)M0 


1.0 


40.0 


3021 


0.4066 


4.51x10"^ 


0.823 


6.20 


1.88 


(1O*,1O*)M0 


1.0 


220.0 


29569 


0.3852 


3.54 xlO"" 


0.924 


6.96 


1.68 



for the former one is 3 times slower than the latter. Since the -IPN correction term 
in the phase is proportional to v~'^, this contribution is about 10 times larger for the 
former case. These 2 contributions make the former constraint greater by more than 
1 order of magnitude. Next, we compare the constraint on wbd which is calculated 
from the inverse of the Fisher matrices. Prom the table, we understand that the 
correlation between other parameters for the former case is weaker by more than 1 
order of magnitude. We think that this comes from the difference in the width of 
effective frequency range of observation. From the table, we see that this frequency 
range is larger for the former case by more than 1 order of magnitude, which solves 
the degenaracies between parameters. 

In Table IIIH we show the pattern-averaged results for massive graviton the- 
ory with (10^ 10^)Mq, (10^ 1O5)M0, (10^ W^)Mq, (IQS, W'^)Mq and (10^, 10'^)Mq 
BH/BH binaries. In this case, larger mass binaries give stronger constraints on Ag. 
This is explained from the correction to the phase velocity Vph which is given asr^ 

Since more massive binaries generate lower frequency gravitational waves, these bi- 
naries give larger corrections and make the constraints stronger. However, if we 
increase the mass parameter too much, the SNR starts to decrease since the dy- 
namical frequency range starts to get narrower, which makes the constraints weaker. 
This is why the constraint on from a (10^, 10®)Mq binary is weaker compared to 
the one from a (10^, W^)Mq binary. In Ref. [IBD, we estimated the constraint from a 
(10'', 1O^)M0 binary with LISA to be 4.06 x 10^° cm. In this case LISA gives stronger 
constraint on Xg than DECIGO. This is because LISA is able to detect lower fre- 
quency gravitational waves than DECIGO where the correction on Vph is greater. If 
we use the extended version of DECIGO shown as thick dashed curve in Fig. [H we 
found that (10^, 1O*^)M0 binary gives a constraint of 4.05 x 10 cm which coincides 
with the one obtained with LISA. This result is obvious since the frequency range 
of this binary is / = 2.36 x 10~^ — 4.00 x 10~^Hz in our case and the noise curves of 
extended DECIGO and LISA are almost identical within this frequency range. 

In Table \TV\ we show the results of error estimation in Brans-Dicke theory for 
(1.4, 10)Mq NS/BH binaries with both pattern-averaged analysis and Monte Carlo 
simulations. The first row shows the ones with pattern-averaged estimate where 
we used only one detector and fixed the SNR to 10. The second and the third rows 
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Table IV. The results of error estimation in Brans-Dicke theory for (1.4, 10)Mq NS/BH binaries. 
The first line shows the results of pattern-averaged estimate. We used only one detector and the 
SNR is fixed to 10. The second and the third lines show the results of Monte Carlo simulations. 
We used two detectors for the analyses and we set SNR= y/200 (corresponding to SNR=10 for 
each detector). We distribute 10^ binaries, calculate the determination error of each parameter 
for each binary and take the average. The second line shows the ones without taking spin 
precession into account, whilst the third line represents the ones including precession, a is 
included in the binary parameters for all the cases. 



cases 


wbd 


A\nM 


Aln-q 


AI3 


Ziln Dl 


Af2s 


Aa 




(10«) 


(10-5) 


(%) 






(10"^str) 




pattern-averaged 


1.342 


0.978 


2.78 


0.190 


0.100 




2.18 


no precession 


0.9774 


1.22 


3.06 


0.186 


1.24 


3.27 


2.15 


including precession 


2.317 


0.350 


0.295 


0.0551 


0.183 


2.52 


0.627 



0.2 



0.15 

'E 

I 0.1 
a. 

0.05 



Fig. 2. The histograms showing the probability distribution of the lower bound of (jJbd obtained 
from our Monte Carlo simulations of 10'' NS/BH binaries in Brans-Dicke theory. We take the 
masses of the binaries as (1.4, 10) M© with SNR= \/200. The thick dotted one represents the 
estimate without precession and the thick solid one shows the one including precession using 
DECIGO. For comparison, we also show the results obtained irP* for (1.4, 1000) M© NS/BH 
binaries with SNR= \/200 using LISA. The thin dotted one shows the one without precession and 
the thin solid one represents the one including precession. The dashed line at ujbd = 2.38 x 10'' 
represents the constraint from PSR J1141-6545.'^ 

represent the results of Monte Carlo simulations with SNR fixed to \/200. The second 
row shows the ones without taking spin precession into account whilst the third row 
shows the ones including precession. We see that inclusion of precession improves 
the constraint on cjbd by a factor of two. In general, binary parameters including 
wbd are correlated with spin parameters /3 and a. When we include precession, we 
obtain additional information about spin which solves degeneracies and reduces the 
estimation errors, making the constraint stronger. Figure [2] represents the histograms 
showing the number fraction of binaries which give the constraint on the Brans-Dicke 
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parameter within each bin of a;BD- The thick dotted one represents the estimate 
without precession and the thick solid one shows the one including precession using 
DECIGO. For comparison, we also show the results obtained iiP^ for (1.4, 1OOO)M0 
NS/BH binaries with SNR=-v/200 using LISA. The thin dotted one shows the one 
without precession and the thin solid one represents the one including precession. 
The dashed line at wbd = 1-47 x 10^ represents the constraint from PSR J1141- 
6545.® We can see that DECIGO can put 300 times stronger constraint than LISA. 
The reasons are the same as the pattern-averaged analysis. 

Unlike the case of LISA, these binaries are thought to be the definite sources for 
DECIGO. The event rate of NS/NS binary mergers is estimated to be 10^ yr~ij3il 
and the rate of NS/BH mergers will be about one order of magnitude smaller than 
that of NS/NS mergers (see Shibata et aZ.^^' and references therein). Therefore it is 
possible to put even stronger constraint by performing a statistical analysis. Defining 
the variance of parameter a) = from i-th binary as cTj, the total variance a is 
given by 



where AT = 1 yr represents the observation time, oq represents the current scale 
factor, r{z) is the comoving distance to the source, h{z) shows the NS/BH merger 
rate at redshift z and r is the proper look back time of the source. aQr(z) and 
g are given ao r{z) = H^^ dz' j ^ n^{\ ^ z'Y ^ Qa, dT/dz = {Hq{1 + 

z)\J /^m(l + ^'Y" + ^a\^^ ^ respectively. Following Ref . [M|) . we adopt the model for 
n(z) given by n(z) = hQR(z), where no = 10"^ Mpc~^ yr~"^ is the estimated merger 
rate today and R{z) = l + 2z (for z < 1), (3/4)(5 - z) (for 1 < z < 5), (for z > 5) 
encodes the time-evolution of this rate. This model gives the merger rate of 10^ yr~^. 
For simplicity, we assume that all NS/BH binaries have the same typical masses of 
(1.4, 1O)M0. We first calculate the variance a{z) for various z using pattern-averaged 
estimate and obtain the total variance a from Eq. The Fourier component of 
the pattern-averaged waveform is given by Eq. (29) of Ref. llSp . To take into account 
the effects of redshift, all we have to do is to replace the masses with the redshifted 
ones: rriNs — ^ (1 + z)mf^s and meH — ^ (1 + z)mBK- The distance in this expression 
is to be understood as the luminosity distance given by Dl = ao(l + z)r{z). 

From the pattern-averaged analysis, we find that observation of 10^ binaries 
with (1.4, 1O)M0 can put a constraint wbd > 2.18 x 10^. This is 94 times stronger 
than the one from a single binary placed at 17 Gpc (corresponding to a binary of 
SNR=10). We calibrate the result of this analysis by using the results of Monte 
Carlo simulation to yield ujbd > 3.77 x 10^. This is 4 orders of magnitude stronger 
than the current strongest bound. 

In Table |Vl we show the results of error estimation in massive graviton theories 
for (10^, 1O^)M0 BH/BH binaries at 3Gpc with both pattern -averaged analysis and 
Monte Carlo simulations.. Again, we have chosen an optimised mass parameter 
which can be understood from Table UTTl The meaning of each row is the same as in 
Table IIVI We can see that the constraint on increases by a factor of two when we 
include precession. Figure [3] represents the histograms showing the number fraction 
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Table V. The results of error estimation in massive graviton theories for (10®, 1O^)M0 BH/BH 
binaries at 3Gpc. The meaning of each line is the same as in Table HVl 



cases 


SNR 




AlnM 


Zilnr; 


AI3 


Aa 


AOs 






(lO^^cm) 


(%) 








(str) 


pattern-averaged 


2044 


1.270 


1.19 


1.69 


9.10 


2.44 




no precession 


2601 


1.266 


1.16 


1.64 


8.92 


2.40 


1.16 


including precession 


2666 


3.349 


0.314 


0.0388 


0.0612 


0.529 


0.0248 




Fig. 3. The histograms showing the probability distribution of the lower bound of Xg obtained 
from our Monte Carlo simulations of 10* BH/BH binaries in massive graviton theories. We take 
the masses of the binaries as (10®, 1O^)M0 at 3Gpc. For comparison, we also show the results 
obtained irPI for (10^ 1O®)M0 BH/BH binaries using LISA. The meaning of each histogram is 
the same as in Fig. [S] 



of binaries that give the constraint of each Xg. For comparison, we also show the 
results obtained iiP^ for (lO"^, 10^)Mq BH/BH binaries at 3 Gpc using LISA. SNRs 
for the gravitational waves from these binaries correspond to 1600. The meaning of 
each histogram is the same as in Fig. [2J We can see that the effect of precession is 
larger for LISA. We expect that this is because the effective frequency range is wider 
for (10^, 1O^)M0 BH/BH binaries with LISA. The lower bound on Xg obtained by 
DECIGO, 3.35 x lO^^cm, is one order of magnitude smaller than that obtained by 
LISA. However, this is still three orders of magnitude larger than the one obtained 
by the solar system experiment. These results show how powerful DECIGO can be 
in constraining alternative theories of gravity. 

Recently, Arun and WilP^ have estimated the constraint on Xg including higher 
harmonics in the waveform. Since they do not include spins, it is important to 
include both higher harmonics and spin precession, and estimate the constraint on 
various alternative theories of gravity. This is left for future work. 
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